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Abstract 

A permutationally invariant n-bit code for quantum error correction can 
be realized as a subspace stabilized by the non-Abelian group Sn- The code 
is spanned by bases for the trivial representation, and all other irreducible 
representations, both those of higher dimension and orthogonal bases for 
the trivial representation, are available for error correction. 

A number of new (non-additive) binary codes are obtained, including 
two new 7-bit codes and a large family of new 9-bit codes. It is shown 
that the degeneracy arising from permutational symmetry facilitates the 
correction of certain types of two-bit errors. The correction of two-bit errors 
of the same type is considered in detail, but is shown not to be compatible 
with single-bit error correction using 9-bit codes. 
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1 Introduction 



Quantum error correction is now well-developed in the case of those stabilizer 
codes jU Ej, which arise as subspaces stabilized by Abelian subgroups of the 
Pauli group. These codes, also known as "additive codes," can be regarded as 
an extension of classical binary codes over to codes over GF(4) which satisfy 
an additional orthogonality condition. They generalize the classical notion of 
distance and thus seem best suited to situations in which all one-bit errors are 
equally likely and the noise is uncorrelated. 

There are other approaches to fault tolerant computation which use structures 
which are resistant to decoherence, e.g., topological quantum computation and 
decoherence free (DF) subspaces or subsystems. (See jTT] and jTTj respectively for 
further discussion and references.) Some physical implementations may also be 
designed to protect against certain types of errors. Much of the current analysis 
is based on simple models using independent errors. In more realistic models 
some types of correlated errors may be more probable than arbitrary two-bit 
errors (and possibly even than certain one-bit errors). Hybrid approaches to 
fault tolerance which combine resiliency (either through encoding or hardware 
design) with error correction may require codes with properties different from 
those stabilized by Abelian subgroups of the Pauli group. 

It is now known [201 that other types of quantum codes, often called "non- 
additive," exist. Although some attempts [21] have been made to develop classes 
of non-additive codes, much of this work, e.g., [T^ I12j. has been for non-binary 
codes. In this paper we consider a natural generalization of stabilizer codes to bi- 
nary codes associated with the action of non- Abelian groups. We concentrate our 
attention on the symmetric group as a case study, and call a code on which the 
symmetric group acts trivially permutationally invariant. We will be particularly 
interested in the use of higher dimensional representations for the correction of 
two-bit errors, and the ways in which the degeneracy associated with permuta- 
tional invariance of code words allows the correction of more two-bit errors than 
would be expected by simple dimensional arguments. 

We find a number of new codes. In particular, we give two new 7-bit codes 
which are impervious to exchange and can correct all one-bit errors together with 
some rather special two-bit errors. We show that the classical 5-bit repetition 
code can correct more two-bit quantum errors than those associated with a single 
type of one-bit error. We show that there is a large family of permutationally 
invariant 9-bit codes in addition to the simple one found in [22] • Unfortunately, 
none of these 9-bit codes is as powerful for two-bit error correction as one might 
expect. 

Although the discovery of new codes is always of interest, we emphasize that 
our primary goal is to study permutationally invariant codes as examples of codes 
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obtained from the action of a non-Abelian group. These non-Abelian groups will, 
typically, be more general than subgroups of the Pauli group. 

It is worth pointing out some significant differences between our approach and 
the "Clifford codes" associated with "nice error bases" as proposed by Knill jjll] 
and developed by Klappenecker and Rotteler [T21 US)- Their approach, which 
considers generalizations of the Pauli group for d > 2, yields non-stabilizer codes 
only for c? > 4; we obtain new non-stabilizer codes for d = 2. (Although our 
approach could, in principle, be applied for any d, we study only d = 2.) In the 
KKR approach, the code is associated with a normal subgroup of an error 
group, but need not come from bases for the trivial representation of A^. We 
retain the requirement that a code subspace is spanned by bases for the trivial 
representation of a group, but the non-Abelian group defining our code need not 
be associated with an error group in the sense of Knill 14j. From a formal point 
of view, our group and error sets reside in an operator algebra associated with 
the usual Pauli group, but we do not use this structure. 

It was recognized earlier PQCHl, in the context of DF subspaces, that quantum 
error correcting codes can be obtained as stabilizers for non-Abelian groups. How- 
ever, the use of higher dimensional irreducible representations for error correction 
was not explored. Moreover, the original philosophy underlying the DF approach 
to fault tolerant quantum computation, namely, to avoid anything which might 
perturb the system out of the stable subspace, is antithetical to active error cor- 
rection. In [21 E] the use of encoding to facilitate universal computation, rather 
than error correction, was introduced. Another important development was the 
generalization of DF subspaces to DF subsystems |jl6j, in which the code itself 
can transform as a higher dimensional representation. The notion of stabilizer 
was then modified in ^21 121] to encompass DF subsystems as well. There is now 
an extensive literature on various aspects of both DF subspaces and systems, in- 
cluding proposals for hybridization of DF methods with active error correction, 
and scenarios in which DF encoding can replace active error correction. We refer 
the reader to J7j for references and further discussion. 

Although motivated by the expected utility of codes capable of correcting 
specific set of correlated errors, we do not present a physical model leading to 
such sets. We deal only with construction of codes, leaving their application 
within a full-fledged scheme for fault tolerance for further investigation. 

The rest of this paper is organized as follows. In the next section, we outline 
the basic set-up and notation we will use. We describe different classes of condi- 
tions associated with one-bit errors in Section 3 and analyze them in Section 4. 
In Section 5, we consider two-bit error correction. In Section 6 we first consider 
some explicit examples of codes for n = 5, 7 or 9; we then show that none of the 
9-bit codes can correct all double errors of one type. 
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2 Preliminaries 



2.1 Stabilizers and error sets 

In the general situation, we have a set of errors £ = {ei, 62 ■ ■ ■ cm} which we want 
to correct. We will also have a unitary group Q which acts on the vector space 
C^". The elements of both S and Q will be linear operators which act on C^". 
Typically, these will be non-trivial linear combinations of elements of the Pauli 
group, rather than simply tensor products of Pauli matrices. In particular, we 
can consider Sn as the group generated by the exchange operators E^s which can 
be written as 



where X^, Y^, Zf. denote the action of the a^, Cy and operators on bit k. Note 
that the set {Eis : s = 2 . . .n} suffices to generate the group Sn- 

Since C^" is invariant under the action of Q, it can be decomposed into invari- 
ant subspaces corresponding to irreducible representations of ^. As is well-known 
|H1 EH! , those subspaces corresponding to inequivalent representations are orthog- 
onal, and those for equivalent representations can be chosen orthogonal. We want 
to exploit the freedom in the latter to construct codes with particular properties, 
and use the additional orthogonality from inequivalent representations for error 
correction. 

Now suppose there is a subspace T which is stabilized by Q, in the sense 
g\w) = \w) for all g E Q and all \w) G T. Consider a subset of errors S' which 
is invariant under Q in the sense gepg~^ G £' for Cp G £'. Then the space S'(T) 
spanned by {ep\w) : Cp G S', \w) G T} is also invariant under Q since 



Hence the space S'{T) can be decomposed into an orthogonal sum corresponding 
to irreducible representations of Q. Since the span of S' itself is invariant under 
Q, it too can be decomposed into a sum of irreducible subspaces. In fact, one can 
regard the two spaces S'{T) and span{£^'} as being decomposed in parallel into 
orthogonal sums corresponding to irreducible representations of Q. 

For example, the set of single bit flips S'x = {Xi, X2, ■ ■ ■ Xn} is invariant 
under Sn- In fact, its span is isomorphic to the standard n-dimensional rep- 
resentation of Sn, which decomposes into the sum of the trivial representation, 
spanned by J2r -^r, and the (n— l)-dimensional irreducible representation, spanned 
by Xi — X2, - - - ,Xi — Xn- Similar considerations hold for the errors {Yi, . . . 
and {Zi, . . . Zn} 




(1) 




(2) 
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The resulting linear combinations of errors Xp — Xq may not be invertible. 
However, this poses no problems for error correction because we will only need to 
"invert" when a measurement shows we are in [Xp — Xg)T which is orthogonal 
to the null space of {Xp — Xq). 

In view of their role as bases for the trivial representation, it is useful to define 
the average errors X, F, Z as 

x = -Vx, F=-Vn z = -y^z, (3) 

n ^ — ^ n ^ — ^ n ^ — ^ 

k=l k=l k=l 

Note that Xr—Xs = (Xi—Xs) — {Xi—Xr), and recall that a code that can correct 
errors in a set S can also correct any complex linear combination of these errors. 
Thus, the error sets 

S = {I,Xk,Yk,Zk,k = 1 . . .n}, and 

S = {I,X,Y,Z,Xi-Xk,Y,-Yk,Zi-Zk,k = 2...n} (4) 
are equivalent. 



2.2 Notation 

The 2"dimensional complex vector space C^" has an orthonormal basis {\v) = 
\vi,V2 . . -Vn)]} indexed by binary n-tuples v = (fi,f2 . . . f„) G (Z2)". If an or- 
thonormal basis |0), |1) for is fixed, this is simply the basis of tensor products 
of the form |f 1) ® |f 2) ® • • • ® |f n) with each t>j G Z2. The symmetric group Sn acts 
on C^" via a natural action on these basis vectors; if V takes [1 . . .n) \—>- (zi . . . z„), 
then V\vi,V2 ...Vn) = |fn,fi2 • • • Wi„)- Define 

Wfc = span{ |t>) : wt(f) = k} (5) 

where wt(f ) is the number of k for which Vk = 1- (This is the classical Hamming 
weight of V.) Then C^" = 0^=0^^ orthogonal direct sum of the Wfc. 

Moreover, each Wfc is invariant under Sn and can be further decomposed into 
an orthogonal sum of spaces affording inequivalent irreducible representations of 
Sn- This yields an orthogonal decomposition of C^" into irreducible subspaces. 
However, unlike the regular representation, some irreducible representations oc- 
cur more than once in C^", and others not at all. Appendix lEl describes the 
decomposition of into irreducible subspaces for n = 5, 7, 9. 

Each Wfc contains the trivial representation, for which we introduce the basis 
vector 

Wk= Yl l^) = 5Z^li^^:^) (6) 

•wt{v)=k V K n—K 
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where the second sum ranges over those permutations V which yield distinct 
vectors \v). Thus (M^, ^k) = {^) ■ Occasionally we will use the normalized vectors 

Wk = (^) Wk- Although normalized vectors are useful for many purpose, those 
denoted Wk are more convenient in combinatoric computations. 
Finally, we will make repeated use of the combinatoric identity 

fN\ f N \ N-2K + J fN+J\ 
[k)-[k-j)= N+J [ K ) 

which holds for J = 1,2 and is easy to verify. We will occasionally use the con- 
vention that (^) = when N < K. 



2.3 Codes 

Given a (possibly non-Abelian) group we define a code C as a subspace of C^" 
which is stabilized by Q in the sense 

g\v) = \v) for all g eQ and all \v) e C. (8) 

If C has dimension 2™, then one can effectively encode m logical binary units in 
n physical qubits. We will restrict ourselves here to the simple case of 1 to n 
encoding, for which m — 1 and C is two-dimensional. A code is often specified 
by an orthonormal basis for C, in which case each basis vector, or "code word" 
can be regarded as a basis for the trivial representation of ^. In the case of two- 
dimensional codes, we can interpret these basis vectors as a logical and 1, and 
will label them \cq) and |ci) accordingly. 

We now consider two-dimensional codes for the group Q = Sn- li \v) — 
\vi,V2 ■ ■ - Vn) is a basis vector of C^" of weight k, (or, equivalently, a binary n- 
tuple of weight k) then {g\v) : g G Sn} is the set of all basis vectors of weight 
k. Therefore, any vector satisfying g\v) = \v) for all g G Sn must have the form 

ttkWk, so that we can write a permutationally invariant code as a pair of basis 
vectors of the form 

|co) = ^afeW^fe and \ci) ^J^^^Wk (9) 

k k 

for some complex numbers and bk with '^,^akbk = 0. A vector X^fe "^fc^fe 
called even (resp. odd) if dk is nonzero only when k is even (resp. odd). 

Note that we have defined a code so that the individual basis vectors are 
permutationally invariant. This is a stronger requirement than that the subspace 
defined by the code is invariant under Sn- However, the distinction is unlikely 
to matter in practice. In the case of two-dimensional codes, the two types of 
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invariance are equivalent whenever n > 3. In general one can have an invariant 
subspace of dimension 2"^ only if it can be written as a direct sum of irreducible 
subspaces whose dimensions sum to 2*"; this will usually consist of 2™ copies of 
the trivial representation, in which case the code words are also invariant. 

We will be primarily interested in codes of the form Q which also satisfy the 
following two conditions (which together imply that n is odd). 

I) bk = an-k or, equivalently, |ci) = (®jXj)|co). 

II) Co is even and Ci is odd or, equivalently, {®jZj)\ce) = (— 1)^|q). 

When (I) and (II) both hold, we can write 

(n-l)/2 (n-l)/2 

|Co) = ^ a2jW2j, |Ci) = ^ an-2j-lW^2j+l- (10) 

j=0 j=0 

In addition to simplifying the analysis and ensuring that certain inner products 
are zero, these assumptions serve another purpose. They ensure that the logical 
X and Z operations can be implemented on the code words by ^jXj and ^jZj 
respectively. Since the actual use of codes in fault tolerant computation requires 
a mechanism for implementing gates on the code words this is an important 
consideration. Moreover, there is little loss of generality in this assumption. The 
operators ®jXj and ®jZj lie in the commutant of Sn- Therefore, they necessarily 
map invariant subspaces of Sn to invariant subspaces of Sn- In the case of the 
code space, we require the stronger condition that 

[®jXj]C = C and [0jZj]C = C. (11) 

When fill!) holds, there is no loss of generality in assuming (I) and (II). These 
simply restrict the choice of basis in way that is convenient and can always be 
satisfied. 

Our goal is to construct a permutationally invariant 2-dimensional code that 
can correct all single qubit errors, and to examine the types of two-bit errors that 
can be corrected. 

As noted at the end of Section d non-Abelian stabilizers were considered 
previously in the context of DF subspace codes. Conversely, one can consider a 
permutationally invariant code as a DF subspace which arises from the highly 
idealized situation in which a quantum computer is completely insulated from 
its environment, but the qubits are the spin components of identical particles 
which interact.^ Then, as discussed in [221, even in the absence of spin-spin 

^For a complete analogy to DF subspaces, only those particles within each logical encoded 
unit would be permitted to interact with each other. Interactions between particles in different 
units, would lead to exchange errors between units. However, since these would appear as single 
bit errors in each logical unit, they would also be correctable. 
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interactions, the Pauli principle induces an effective interaction between qubits 
whose DF subspace group is precisely that generated by exchanges. In essence, 
the Pauli principle requires correlations between the spatial and spin components 
so that spatial interactions (such as the Coulomb interaction) affect the spin 
components. The result of tracing over the spatial component yields a completely 
positive map on the spin components, as in the standard noise model. Only 
a fully symmetric spin function allows the full wave function to be a product 
(with an anti-symmetric spatial function) consistent with the Pauli principle for 
fermions. Thus, a DF subspace is precisely one which transforms as the trivial 
(or fully symmetric) representation of Sn- Although this is not a very realistic 
DF scenario, it is useful to see how codes constructed for different purposes can 
be interpreted within the DF subspace , as well as the stabilizer, formalism. 

3 Error correction conditions 

The now well-known necessary and sufficient condition [3 for the code C to 
correct errors in a set £^ = {ei . . . cm} can be stated as 



where the matrix dpq does not depend on One often chooses codes for which 
dpq = Spqfip, but that is not necessary. Indeed, the requirement 



which is implicit in fll2|) . implies that one can always transform the error set into 
a modified one S for which the stronger condition dpq = 6pqfip holds. 

Strictly speaking one can only determine whether or not a particular set of er- 
rors is correctable; not whether a particular error or type of error is "correctable" . 
However, it is often natural to look for codes for which the set of correctable errors 
includes all errors of a particular type, e.g., the one-bit errors. One can then ask 
what additional errors could be added to this subset to yield a set S satisfying 
(fT^ . In our discussion of such situations, the subset involved may be implied by 
the context. 

We will find it useful to think of (|12j) as defining a pair of matrices D** with 
elements d^^ = {epCi,eqCi) for 2 = 0,1 and a matrix B = D^^ with elements 
bpq = dp^ = {epCo,eqCi). Then flT^ is equivalent to the requirements B = and 
Z^oo = 1)11 (Because d^^ and d^^ are complex conjugates, D^^ = <^ D^'^ = 0. 
Hence we need not consider D^^ explicitly and will use only B = D^^.) For 
simplicity, we omit the superscript in dl . 





(12) 



dpq = {epCo,eqCo) = {epCi,eqCi) 



(13) 
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For example, the 9-bit permutationally invariant code in j22l corrects single 
qubit errors as well as the Pauli exchange errors (transpositions) E^s, for the 36 
unordered pairs r, s. We can consider the above matrices with respect to the errors 
S = {/, Ers, Xi, . . . , Xg, Fi, . . . , Fg, Zi, . . . , Zg} . It was shown that = D^^ and 
has the block diagonal form 



(14) 



where is a 37 x 37 rank one matrix and the 9x9 matrices Dxx, Dyy, Dzz 
correspond to the one-bit errors indicated by the subscripts. These all have the 
cyclic form 



[Do 








\ 





Dxx 














Dyy 













DzzJ 



(a b 
b a 

\b b 



b 

a J 



(15) 



For any permutationally invariant code the blocks Dxx, Dyy, Dzz necessar- 
ily have the form (fT3j) . Such matrices can always be diagonalized by a change 
of basis to (1, 1 ... 1) and its orthogonal complement. This corresponds to re- 
placing the errors {/i, /2 . . . /«} by the corresponding average {/} and a suitable 
orthogonalization of {fi — fk, k = 2 . . . n} where / denotes any of X,Y, Z. 

Now the orthogonality of subspaces associated with different irreducible rep- 
resentations ensures that 



{fW„{gr-gs)Wk) = 



(16) 



for all j, k and any choice of / = /, X, Y, Z and g = X,Y, Z. Alternatively, we 
can show this directly by observing that the exchange operator Ers is unitary so 
that 

{JWj, {gr - gs)Wk) = {Er.jErs{ErsWj),Ers{gr ' g s) Ers{ErsWk)) 

= (fW^A9s-gr)Wk) 
= -{fW„{gr-gs)Wk) 



which implies ()16p . For such codes, each of the matrices D^^, {i = 0, 1) and B 
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have the form below (which we write only for D) with respect to the order in P)). 



/ dii 


dix 


djY 


diz 


dxi 


dxx 


dxY 


dxz 


dyi 


dyx 


dyy 


dyz 


dzi 


dzx 


dzY 


dzz 











V 



Dxx 


DxY 


Dxz 




Dyx 


Dyy 


Dyz 




Dzx 


DzY 


Dzz 


J 



(17) 



Conditions (I) and (II) immediately give many additional zero entries. One 
nice way to see which entries are zero is to observe that ^kZk commutes with 
Zr and anti-commutes with Xr and Yr for all r. Thus, for every one-bit error Cp, 



zZk) 



:Zk)ep, where 

+1 for ep e {lX_{Zr - Z,)} 

-1 for epe{X,Y,{Xr-X,),{Yr-Ys)} 



Also, 



)kZk is unitary so that 



{{®kZk) epCi, {®kZk) egCj) 
= e^eg {ep{®kZk)ci,eq{^kZk)cj) 

From this we can conclude the following. 

A) When i = j, {epCi,eqCi) = whenever ep ^ Thus, dix = djy = dxz = 
dyz = and Dxz = Dyz = 0. 

B) When i ^ j, {CpCi, CpCj) = whenever = e^, from which we can conclude 
biz = bzi = hxY = byx = and Bxx = Byy = Bzz = Bxy = Byx = 0. 



Combining this with djg 
that the have the form 



^ ^ dgf = ^ and Dfg = -^^ Dgf = 0, we find 



/ dn 








diz 







dxx 


dxY 










dyx 


dyy 







dzi 








dzz 












Dxx 











Dyx 














DxY 
Dyy 







Dzz ) 



(19) 
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and B has the form 






bix 


bjY 





bxi 








bxz 


byi 








byz 





bzx 


bzY 














V 



B 









Bxz 








Byz 


zx 


BzY 


) 



(20) 



Now, observe that ®kXk commutes with and anti-commutes with and 
Zr- Proceeding as above, we find 



I for epG{/,X,(X,,-X,)} 

4 for {F,z,(i;-y;),(z,-z,)} 



(21) 



and 



^ ef (ep (OfcXfc) Q, eg (®fcXfc) 9) 



XX/ \ 



where we interpret i + 1 and j + 1 mod 2. Thus we can conclude 

C) When i = j, condition (fT^ holds whenever = e^. (This means, in 
particular, that the diagonal entries and blocks of D^^ and D^^ agree.) 

D) When i = j and 7^ e^, condition (|T^ can only be satisfied if (cpQ, CgCj) = 
for z = 0, 1. Thus we must have djz = dxy = and Dxy = 0. 



E) When i 7^ j, (epCo,egCi) = ±(epCi,egCo) = ±(eqCo, CpCi). Thus, we can 
conclude, e.g., that matrix entries bxz = -v^ bzx = and blocks Bxz = 
<S=^ Bzx = 0, so it suffices to check entries of B above the main diagonal. 

Thus, when conditions (I) and (II) are satisfied, we find that sufficient (and 
necessary) conditions for (fT^ to hold are that 

• All off-diagonal entries and blocks in (|T9|) are zero, 

• All remaining entries in ()2()|1 are zero. 

Moreover, it suffices to check matrix elements above the main diagonal in ()19|) 
and (1201). 

We can break these conditions into several groups, which will turn out to be 
related or equivalent. 
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a) Conditions bjx = bjY = bzx = bzY = 0, which are equivalent to 
(co, {X ± iY)ci) = {Zcq, {X ± Y)ci) = 0. These will yield just two con- 
ditions when the are real. 

b) Conditions djz = dxv = 0, i.e., {cq.Zcq) = {Xcq^Ycq) = 0. These will 
reduce to one condition when the are real. 

c) The block conditions Bxz = Byz = 0, which are equivalent to 
{[{Xi - Xr) ± i{Yi - Yr)]cQ, (Zi - Zs)ci) = for 2 < r, s < n. 

d) The block condition Dxy = 0, which is equivalent to 

((Xi - Xr)co, iYi - r,)co) = for 2 < r, s < n. 

We will see that for codes which satisfy conditions (1) and (11) and have all 
coefficients real, conditions (c) on blocks will be satisfied whenever (a) holds; and 
conditions (d) on blocks will be satisfied whenever (b) holds. Thus, we will only 
need to satisfy three non-linear equations for such codes. We can summarize this 
as follows. 

Theorem 1 Assume that the coefficients ak associated with a permutationally 
invariant code which has the form lll(J\} and length n are all real. Such a code can 
correct all one-bit errors if and only if the following equations hold. 

n + 1 I n 



( ri+l ) "ll±l +2 ^ a2man-2m+l ( 2^ ) 
\ 2 / ^ m=l ^ ^ 

f n\ L(«-3)/4j , ^ 

( ] aizJ. +2 ^ a2man-2m-i {u - 2m) ( 

V 2 / rr>— n V 



2 / m=0 

5Z ''2m(«-4m)( 

m=0 ^ ^ 



(23) 
(24) 



The theorem will follow from the analysis in the next section. The result can 
be extended to complex ak as discussed in Appendix O As noted before, it is 
implicit in (|Tn|l that n is odd. 
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4 Error condition analysis 

4.1 Conditions of type (a) — co,ci orthogonality 

First we give expressions for the action of the average errors X, Y and Z. It is 
easy to check that for < k < n 

nZWk = {n-2k)Wk (25) 
nXWu = {k+l)Wk+i + {n~k+l)Wk-i (26) 
inYWu = {k + l)Wk+i - {n - k + l)Wk-i, (27) 

where it is understood that if m < or m > n, Wm should be replaced by zero; 
for example, nX Wq = Wi. 

To analyze the requirement bix = biy = bzx = bzv = 0, it is equivalent (and 
somewhat easier) to use the conditions 

= (co, (X±zF)ci), and (28) 
= (Zco,(X±iY)ci) (29) 

When condition (I) holds, these conditions become 

n(co, (X + ^Y)c^) = a^a^-k+i = (30) 

k=i ^ ^ 

n^Zco, (X + iY)ci) = 2 - '^k)k Q aka^-k+i = (31) 

n(co, (X - tY)ci) = 2 ^ (n - A;) ( ^ j Ukan-k-i = (32) 

k=0 ^ ^ 

n-1 y N 

n^Zco, (X - tY)ci) = 2j2in-'^k){n-k)r\akan-k-i =0 (33) 

k=0 ^ ^ 

Thus far, condition (I) has played a minor role and one can easily obtain more 
general conditions by replacing a„_fc by bk above. Now, however, we make explicit 
use of the fact that all products have the form a^an-k+i to conclude that the real 
parts of the expressions in (jHUI) and and in and agree up to sign, 
which leads to the following lemma. 

Lemma 2 When the coefficients are all real, the equations / t.vrjj) to are 
equivalent in pairs, UIS) ^ /f371 ) and ^EE)- When n is odd these reduce 
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to 

(n-l)/2 



QfcQn-fc+l2fc ( 7 j + '^"+1 ( n+l 

fe=l \ / ^ 

afca„_fc_i2(n - /c) I , ) + 



2 n+l 

2 I n-1 



(34) 



(35) 



When n is even, similar expressions hold with upper limits of n/2 and n/2 — 1 
respectively hut without any extra square terms analogous to a\±i . 



When condition (II) holds, equations (jMj) and (jH^j) reduce to and 
which proves the first part of Theorem ^ 

Proof: First, observe that for m 7^ the term aman-m+i occurs twice in (jSOI), 
once for m = k and once for m = n — k + 1. Thus, the coefficient of amOn-m+i is 



mi )+(n — m+l)( ) = 2m 

\mj \n—m + l 

The coefficient of the same term in (|31|) is 



\m J 



in — 1m) m 



( ^ I + fr2-2(n-m + l)l (ra-m+1) ( ^ ) = -2m ( 
\m) ^ ^ ' \n-m + lj \ 



To obtain a general proof and reduction to (j22p . it suffices to make the change of 
variable k —>■ n — k + 1 when k > in (jHIH ) and (jSH) and, as above, use the 
elementary identity 

ml^]=(n-m+l)l ^ ]. (36) 
\mj ^ '\n-m+l) ^ ' 

Similarly, the change of variable fc— — fc — lforA;> in (jH^ ) and ()33|) 
yields (j22I). QED 

4.2 Conditions of type (b) — off-diagonal conditions 

Using (j2SI) to (j2I|) one finds that 



l,nXW^,%nYWk) = {W,, nZW^) = (n - 2k) 
{nXWj, i nYWi>) = if j - ^ ^ ±2, 0. 
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From these relations, it follows that when all and hk are real 



n^Xc,, I Yc,) = n{c,, Zc,) = " 2^^) (?) (37) 

k ^ ^ 

where ak equals ak or bk according as j equals or 1. Thus, we can conclude that 
when ak,bk are real, the off-diagonal conditions djz = dxY = hold if and only 
if 



{31 



This reduces to to (|24j) when conditions 1 and (11) are satisfied. If some a^, bk 
are not real, then (jHTj) holds with al replaced by la^p or but the additional 
condition ()92|) is needed to ensure that the imaginary part of {Xci,iYci) is zero, 
as discussed in Appendix O 

4.3 Conditions of type (c) — block co,ci orthogonality. 

It will again be useful to replace the separate X, Y equations by their sums and 
differences. The requirement that the blocks Bxz = Byz = is equivalent to 

( [(Xi - X,,) ± i{Y^ - Yr)] Co, (Zi - Z,) ci) = 0. (39) 

for 2 < r, s < n. We now need results from Appendix A. When conditions (1) 
and (11) hold, equations (j8na|l . (j8nb|l and ()8nc|l imply that (|39|l is equivalent to 
the following pair of equations 

("-3)/2 , _ 

a2ma'n-2m^i\ „ ) = 0) ^ud (40) 

V 2m / 

m=o ^ ^ 

a2ma„_2m+i f „ _ „ ) = 0. (41) 



m=l 



To see that these are equivalent to (j^^ and ()23|1 , again make a change of variable 
of the form A;— >n — fc^lin the second half of each sum and use the identities 

n — 2\ , n — 2 \ fn — 1\ n — k fn 



k J \n — k — 1 J \ k J n \k 

'^^~^''\l-2) + [n-k-l) = [k-l) = Hi 
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4.4 Conditions of type (d) — block ofF-diagonal conditions 

We now consider the condition dxY = which means that 

{{X,-Xr)co,t{Y,-Ys)co) = (42) 

for all choices of 2 < r, s < n. The crucial fact is that the inner products of 
this type with r = s and r ^ s differ only by a factor of 2 as shown by ()78|1 in 
Appendix A. 

Theorem 3 When conditions (I) and (II) are satisfied and the Ok are real, 
holds if and only if l{24\ ) does. 

Proof: It follows from ()79b|) and ()79c|) that when 1, r, s are distinct 
((Xi-X,)co,z(ri-r.)co) = 2((Xi-X,)co,z(Yi-n)co) 



{n-3)/2 



2Ao + 2 f 

m=l ^ 



n-1 
2m 



(43) 



where = |a2mP - |«2m+2p + a2m«2m+2 - a2ma2m+2- Thus when the are real, 



will be satisfied for all choices of r, s if 



(n-3)/2 
2 I \ ^ 2 

m=l 



n - 2 
2m 



- 2 
2(m - 1) 



(n - 2)aLi = 0. 



(44) 



One can then conclude that (ji^ is equivalent to if n ('^^^) — (2(^-1)) 
(ra - 4m) (2^) which follows from ((Tj) with N = n - 2,K = 2m, J = 2. 



5 Two-bit errors 



5.1 Some special types of two-bit errors 

The standard 7-bit CSS code |K| can correct two-bit errors of the form XrZg but 
not those of the form X^Xg or Z^Zg. The last two are far more likely to occur, 
especially for nearest neighbors. We now consider the effect of two-bit errors of 
the same type, which we call "double" errors, on permutationally invariant codes. 

Recall that exchange errors have the form (^Q). Permutationally invariant 
codes are designed so that exchange errors are degenerate with the identity, i.e.. 
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Ers\cj) = |cj) for j = 0, 1. Now consider the following three errors 



rs 


1 
2 












^z. 


(45a) 


Gj-s 


1 

2 






^Xs 


-Yr^ 


^Ys-Zr^ 


^Zs 


(45b) 


Hrs 


1 
2 




^ I - Xr^ 




+ Yr(i 


^Y.-Zrd 


^z. 


(45c) 



and observe that 

• Fj-s exchanges two bits and multiplies by —1 if and only if the values of the 
bits are different. 

• Grs flips the two bits r and s if and only if they are the same. 

• flips the two bits r and s and then multiplies by —1 if and only if they 
arc the same. 

In a product basis of the form |00), |01), |10), |11) these operators are represented 
by the matrices 



/I 








o\ 




(I 








o\ 








1 













-1 








1 








p — 

^ rs — 





-1 
















V 












1/ 


/o 








\\ 




/o 








-1\ 





1 













1 














1 





Hrs 








1 





^1 








0^ 




V-1 








0^ 



Any code which can correct all errors of the type Ers, Frs,Grs, Hrs can also 
correct any error of the form ZrZs, XrXs,YrYs, since an error of one type can 
be written as a linear combination of those of the other. For permutationally 
invariant codes, these two types of errors are actually equivalent. 

Theorem 4 // \ip) is permutationally invariant (i.e., Ers\ip) = \ip) for all r,s), 
then the operators Frs,Grs and Hrs have the same effect on lip) as ZrZg, XrXg, 
and YrYs respectively, i.e., Frs\ip) — ZrZs\ip), Grs\i^) — XrXs]^)) and Hrs\ip) — 

YrYsl^). 

Proof: First note that Ers + Frs — I + ZrZg. Then 

Frs\il^) = {I + ZrZs - Ers) IV') = IV') + ^r^.lV') " = ^r^.lV') 

The other two cases are done similarly using Ers + Grs and Ers + Hrs respectively. 
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5.2 Two-bit error correction conditions 



We begin with some simple, but fundamental, results. The first follows from the 
fact that all double errors preserve parity. 

Theorem 5 Whenever condition (II) is satisfied, 

(cpCo, CgCi) = (46) 

for any pair of errors in the set {I, Z^Z^, XrXg, Yj-Yg} or, equivalently, in the set 

{I , Ersy Frsj Grrs: Hrs\ • 

The next theorem says that all inner products of the form {ZrZgCj, ZgZfCj), 
{ZrZsCj,XqXtCj) etc. are independent of j = 0,1. It follows easily from the 
equivalence of condition (I) to {®kXk) |co) = |ci), and the fact that ®kXk is a 
unitary operator which commutes with any error of the form X^X^, Y^Ys or Z^Zg. 



Theorem 6 Whenever condition (I) is satisfied, 

{frfsCo, QqQt Co) = {frfsCl, ggQi, Ci) (47) 

where /, g denote any of {X, Y, Z} (the same as well as different) and r, s, q, t are 
arbitrary. 

One is often interested in knowing which two-bit errors can be corrected in 
addition to one-bit errors. Conditions involving the average error ZZ can be 
readily calculated by noting that 

J2 ZrZs = J2 ^^^^ ~J2^r = [Yl ^ - = i^'^f - (48) 

r^s r,s r r 

Combining this with one finds 

n{n - l)ZZWk = ( ZrZs) Wk = [{n - 2kf - n] Wk- (49) 

The additional conditions needed to correct all errors of the form Z^Zg include 
(Zco, ZZcq) = and {ZZcq, {X ± iY) ci) = 0. The latter gives the following pair 
of conditions 

= 2j2[{n-2kf-n]k(j^akan-k+i (50) 

n-l ^ N 

= 2 5^[(n-2A;)2-n](n-A;) (^j afca„_fc_i (51) 
fc=0 ^ ^ 
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Using and (jl^ . one finds that {Zcq, ZZcq) = dz,zz = is equivalent to 

= "^\ak\\n-2k)[{n-2kY -n](yj- (^2) 

k ^ ^ 

Although one can write down a similar set of conditions for the correction of 
errors of the form X^Xg, it is probably easier to use the following observation. 

Theorem 7 A permutationally invariant code |co), |ci) which satisfies conditions 
(I) and (II) [and corrects a specified set of one-bit errors] can correct all errors 
of the form ZrZg if and only if the code 

IQ) = //®"|co) + (-l)^'i/®"|ci) (53) 

can correct all errors of the form XrXg. 

The map \cj) t-^ \Cj) in consists of a Hadamard gate acting on all qubits, 
followed by an effective Hadamard operation on the resulting code words them- 
selves. This is extremely useful and is its own inverse. We will refer to it as the 
"Hadamard code map". 

5.3 Degeneracy enhancement of classical codes 

It follows from Theorems El and IHl that every permutationally invariant code for 
which both conditions (I) and (II) are satisfied can correct all double errors of 
the form Zj-Zs, XrXs,YrYs provided that we do not also require single bit errors 
to be correctable. For example, the 3-bit repetition code |co) = |000), |ci) = |111) 
is generally regarded as able to correct all single bit flips, but no other errors. 
However, one could instead use it to correct all double bit flips, at the expense 
of the ability to correct any single bit errors. The theorem above says that it 
can do even more — it can correct all two-bit errors of the same type. Although 
this might seem surprising at first, it is easy to understand why it is true. For 
this code ZrZs\cj) = \cj) so that ZrZg is degenerate with the identity. Similarly, 
YrYg is degenerate with XrXg. Note that this degeneracy extends to any n-bit 
repetition code 

|co) = |00 . . . 0) = Wo, |ci) = 111 . . . 1) = (54) 

When n > 5, the simple repetition code (|54j) can correct all single and all 
double bit flips. Indeed, for n = 5, this is just a classical code for two-bit error 
correction. Applying the Hadamard code map ()53|) to yields a code which 
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can correct all single and double phase errors. In fact, omitting the normalizing 
coefficients, this code is 



\Co) = 




\Ci) = 





m 



(55a) 

2m+l (55b) 



Because the phase errors preserve parity, the necessary and sufficient condi- 
tions for a code satisfying conditions (I) and (II) to correct both single and double 
phase errors are 

= (co,Zco) (56a) 
= (Zco,ZZco) (56b) 
= ((Zi-Zi)co,ZiZ,co) = it^s). (56c) 



Note that since {{Zi — Zt)co, ZZcq) = 0, the single and double-Z errors which 
transform as the (ra — l)-dimensional representation are orthogonal if and only if 
()56c|) holds. In fact, as shown after (jHBj) in Appendix A, ()56c|) is redundant, i.e., 
it is satisfied whenever fl56aj) and ()56b|) hold. Thus, one finds that the necessary 
and sufficient conditions for a code satisfying conditions (I) and (II) to correct 
single and double phase errors are (|HH|l [which becomes (|Hn|) when is complex] 
and (jS21), which we rewrite below. 

= ^|afc|'(n-2A;)('^') = (57a) 
k ^ ^ 

= ^|a,f (n-2A;)[(r2-2A:)'-n](^^^ =0. (57b) 

When n = 5, the pair of equations in (j37|) has exactly one solution (up to 
normalization), namely | flop — l'-^2p — \(^4 

I . Choosing identical phases yields the 
code in (j55|) . In addition to correcting all one and two-bit phase errors, it can 
also correct all errors of the form XrXg and YrYg. Choosing other phases yields 
other codes and taking the Hadamard code map yields classical codes for two-bit 
error correction that are distinct from These also satisfy conditions (I) and 
(II) and, hence, can correct all double ZrZg and YrYs errors as well as single and 
double bit flips when used as quantum codes. 

When n > 7 and odd, the pair of equations (jKTj) has infinitely many solutions 
in addition to |flo|^ = |fl2p. . . = |fln-ip- Taking the Hadamard transform then 
yields infinitely many classical codes for two bit error correction. 
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5.4 Higher dimensional representations 



In this section we take some preliminary steps toward exploiting higher dimen- 
sional irreducible representations for correction of errors in addition to one-bit 
errors. First, we review the mutually orthogonal subspaces required for the cor- 
rection of single errors. The operators J, X, Y, Z acting on the code words |co), |ci) 
require four pairs of one-dimensional subspaces which transform as the trivial rep- 
resentation. The three sets of differences Xi— Xr, Yi — Yr, and Zi — Zj. acting on 
the code words require three pairs of subspaces of dimension n— 1 which transform 
as the even (n— l)-dimensional representation. But (as described in Appendix B) 
the decomposition of C^" into an orthogonal sum of irreducible subspaces includes 
other irreducible representations of Sn- 

The next irreducible representation has dimension and arises in the 

decomposition of double errors of one type, e.g., frs = XrXg. Consider the 
subspace generated by frsWk for A; = 2 (or A; = n — 2) as r, s run through all 
"•^"■'^^ combinations of r < s. This subspace splits into an orthogonal direct sum 
consisting of 

• a 1-dimensional subspace spanned by // Wk where // = (2) ^ J2r=/=s frs is 
the average error of this type, 

• an (n— l)-dimensional subspace spanned by the vectors frW^ for r = 2,3 . . .n 
where = X;"=2 fis - Es^r frs r = 2, 3 . . . n , and 

• an "^"2"^^ -dimensional subspace which is obtained by taking the orthogonal 
complement of the vectors ffWk and frWk in spanj/rsW^} ■ 

As described in SectionEHl the error set {frs} can be correspondingly decomposed 
into bases for representations of Sn with dimensions 1, — 1, and 

For an explicit example of the last type of error, consider n = 4. Then 
VV2 splits into three subspaces, corresponding to irreducible representations of 
dimensions 1, 3 and 2. The last is spanned by the vectors: 

2/12 — /l3 — /l4 — /23 — /24 + 2/34 
fl2 + /l3 — 2/14 — 2/23 + /24 + fsA 

There is a sense in which these errors are rather delocalized, since they act on all 
six pairs of qubits. Although one could eliminate some pairs by a different choice 
of basis vectors, one can not, e.g., eliminate all terms of the form fj^ involving 
the 4th qubit. This delocalization is, unfortunately, the antithesis of what one 
might want in certain situations, such as errors between nearest neighbors. 

We now focus on n = 7 as an example and note that can be decomposed 
into an orthogonal sum of irreducible subspaces spanned by 
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• 8 orthogonal bases for the trivial 1- dimensional representation, 

• 6 orthogonal bases for the even 6-dimensional representation, 

• 4 orthogonal bases for a 14-dimensional representation, 

• 2 orthogonal bases for a second, inequivalent, 14-dimensional representation. 

Thus, correcting the one-bit errors requires all of the available 1 and 6 dimensional 
representations. However, the two types of 14-dimensional representations are 
available to correct two-bit errors. One of these 14-dimensional representations 
has spin | and is associated with the partition [5, 2]. This is the "-("^-3) -dimensional 
representation which arises in the decomposition of frs errors described above, and 
can be used to correct the corresponding subclass of double errors. There are three 
kinds of double errors, those from Z^Zg, from X^X^, and from Yj-Ys. This would 
seem to require six orthogonal subspaces which transform as this 14-dimensional 
representation; however, we have only four — one each from W2,W3,>V4,>V5. 
Nevertheless, Theorems El and El imply that all three types of double errors can 
be corrected. 

This is indeed the case and is the result of degeneracy. For permutationally 
invariant codes, 

{XrX, + ZrZs + YrY^m = \^) . (58) 

Thus, it suffices to correct any two of X^Xg, ZrZs, Y^Yg to ensure that all three 
types of errors can be corrected, and this requires only four 14-dimensional sub- 
spaces, exactly what one has available when n = 7. Thus, a 7-bit permutationally 
invariant code which can correct all one-bit errors can also correct all errors in 
the 14-dimensional irreducible components of the decompositions of XrXg, ZrZg, 
and YrYg. This implies that arbitrary double errors would be corrected about 2/3 
of the time. Similarly, a 9-bit code could correct them about 3/4 of the time. 
Unfortunately, the other 1/3 (or 1/4) of the time, the procedure does not simply 
fail to detect the error — it incorrectly interprets a two-bit error as a one-bit error 
and the attempted correction actually introduces additional errors. 

We next consider the case n = 9. Correcting the one-bit errors uses 8 of the 
10 available 1-dimensional representations and 6 of the 8 available 8-dimensional 
representations. In addition to the six 2 7- dimensional representations, two 1- 
dimensional representations and two 8-dimensional representations are also po- 
tentially available to correct some two-bit errors. Correcting one type of frs double 
errors would require a pair of 1-dimensional, 8-dimensional and 27-dimensional 
subspaces. Based on dimensional considerations, one might expect to correct 
one type of double error completely using a 9-bit permutationally invariant code. 
Unfortunately, as will be shown in Section \b.b\ this is not possible. 
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One can still ask what additional errors can be corrected with permutationally 
invariant 9-bit codes. The operators /, X, Y, Z acting on the code words \cq), |ci) 
generate an 8-dimensional space. Taking the orthogonal complement in the 10- 
dimensional subspace spanned by {Wq, Wi . . . Wg] yields a two-dimensional sub- 
space. There is a family of linear operators which map |co), |ci) to a pair of 
orthogonal vectors in this two-dimensional subspace. Any member of this family 
can be chosen as an additional correctable error. Similarly, there will be a set of 
correctable errors which transform as the 8-dimensional representation and whose 
action on the code words spans the orthogonal complement of the one-bit errors 
in span( Although a procedure for obtaining these operators can be 

written down, we have been unable to characterize them in a useful way. 

6 Special cases 
6.1 n = 5 

When n = 5, conditions (I) and (II) hold, and all at are real, the three necessary 
and sufficient conditions in Theorem ^ become 

02(34 = 

O0O4 + 3o2 = 
Oq + 202 — 304 = 0. 

It is easy to verify that these have no non-trivial solution. This is not surprising. 
It is well-known that the 5-bit code for correcting all one-bit quantum errors is 
essentially unique and is not permutationally symmetric. 

Nevertheless, there is still something to be learned by looking at 5-bit codes. 
As discussed in section 1^751 the simple repetition code 

|co) = lOOOOO), |ci) = llllll) (59) 

corrects both all single and all double bit flips, and 

|Co) = H^'>\cq) + H^'>\ci) = W0 + W2 + W4 (60a) 
\Ci) = H^^\co) - H^^\ci) = W^ + W^ + Wi (60b) 

corrects all single and double phase errors. In fact, when n = 5, equations ()57a|) 
and ()57b|) imply that the only codes satisfying conditions (I) and (II) are those 
with |oop = |o2p = |o4p. 

Moreover, both codes can correct all double errors of the form XjXk, YjYk, 
ZjZk- To see this, note that Z^Zglip) = {ip) on the span of ()59|1 so that implies 
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XrXs\ip) = —YrYs\ip), i.e., the pair {ZrZs,!} is degenerate and this induces a 
degeneracy on the pair {X^-Xg, Yj-Ys}. 

Thus, the 5-bit codes ()59|) and ()60|) can each correct more types of quantum 
errors than one might expect from their classical distance properties. They are 
optimal for the correction of all one-bit and two-bit errors of a particular type 
(phase or bit flip) and can not correct additional one-bit errors. Nevertheless they 
can correct additional types of two-bit errors. 



6.2 n = 7 



When n = 7, conditions (I) and (II) hold, and all are real, the three conditions 
in Theorem n become 



3a2a6 + 5a4 = 
Ooag + 150204 = 



On + 9ao 



0. 



(61a) 
(61b) 
(61c) 



It is not hard to see that oe = implies all = 0. Therefore we can divide 
through by ag or, equivalently, assume without loss of generality that = 1. 



Then fj61a|) and ()61b|) imply 02 
substituting in ()61c|) yields 



5^2 



ai and ao = 25a|. Letting x = ai and 



125x^ + 5x'^ -x-1 



to which X = \ is the only real solution, giving 



ao 



±V5 



1 



0^2 



a4 = ±- 



1 



ag = 1. 



3 75 

It is then straightforward to verify that both signs in the normalized vector 



I Co) 



± VlEWo - V7W2 ± V21Wi + V21W( 



(62) 



yield acceptable codes. This gives two distinct new codes when n = 7. 

It is interesting — and a good check — to write the vectors X\ci), Y\ci), Z\co) 
and see that together with |co) they form an orthogonal set. Up to normalizing 
scalar s, we have 

^|co) = 
X|ci) = 

y\ci) = 



± VsbWo - V3W2 TWi- 51?f 

V7Wo + V3{2 ±V5)W2 + (±4 -V5)W4 + (±^5 - 2) 
-V7Wo + V3(2 ^ VS) t?2 + (±4 + V5)W4 + (tVs - 2) 
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Remark: One might ask if one can obtain additional permutationally invariant 
7-bit codes by allowing complex coefficients. In that case the following equations 



are necessary and sufficient. 

= 10|a4p + 3(02^6 + ^2^6) (63a) 
= 02^6 ~ Ci2(ie (63b) 
= {aoCLe + cto^e) + 15(0204 + 02^4) (63c) 
= (doOe — floOe) + 5(02^4 — 02^4) (63d) 
= |aoP + 9|a2p — 5|a4p — 5|a6p (63e) 

= (aoa2 — 0902) + 10(a2a'4 — 0.20^4) + 5(0405 — 0406)- (63f) 

As in the real case, Og = forces all coefficients to be zero, so we can assume 

Oe = 1. Then (j63b|l implies that 02 is real and (j63b|l that 02 = —(5/3) [04 p. 



However, ()63d|) and ()63fj) yield a pair of linear equations for Imoo and Ima4 
which have a non-zero solution if and only if 02 = +1 which is not consistent with 
02 = —(5/3) 1 04 p. Thus, there are no permutationally invariant 7-bit codes other 
than those in (j62| . 

6.3 n = 9 

When n = 9, conditions (I) and (II) hold, and all ak are real, the three conditions 
in Theorem n become 

a2ag + Ta^a^ = (64a) 
35a4 + aoOg + 280205 = (64b) 
al + 20al + Ual - 28a^ - 7al = 0. (64c) 

We now show that these have have infinitely many solutions. First, suppose 
08 = 0, so that the ffist equation becomes 0405 = 0. If og = 0, then all the 
coefficients are zero. If 04 = then we find the only possibility is 02 = and 
Oq = 28ag, giving the two solutions found in jj22j. To find the remaining solutions 
we may assume ag = 1. 

If Oe = 0, then also 02 = and oq = —3 504 where 04 = x is the positive root 
of the quadratic equation 175x^ + 2x — 1 = 0. This gives two more solutions. If 
ao = 0, then we find all of the remaining coefficients depend on Og = t, where t is 
a positive root of the cubic f{t) = {28^/5)t^ + (2 ■ 28V5)t^ - 4t - 1: 

tti = (28/5)t, as = ±Vi, 02 = (-7 ■ 28)t. 

Since /(O) < and /(I) > 0, / does have a positive root (approximately t = 
0.478), and this gives two more solutions. 
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There are no further solutions with any of the equal to zero, so now assume 
they are all nonzero. Writing x = a\ and t = 04, it follows that 

02 = ±7^/xt ao = 7(28tx - St^) 

where x and t satisfy the equation 

x^{5A88t^) + 4a;(-490t^ + 35t^ - 1) + (175t^ + 2t^ - 1) = 0. 

Using Maple, one can verify that there are infinitely many values of t (e.g., all 
for which —0.25 < t < 0.4) for which this quadratic in x has at least one positive 
solution. Thus there are infinitely many solutions in which all the coefficients 
are non-zero real numbers. 

6.4 Conditions for double error correction with 9-bit codes 

As discussed in Section 15.21 a 9-bit permutationally invariant codes which can 
correct all errors of type Z^Zg as well as all one-bit errors, must satisfy at least 
9 conditions. In the notation of Section El there are six of the form bjx = bjy = 
bzx = bzY = bzz,x = bzz,Y = and three of the form dxY = djz = dz,zz = 0. 

First, consider the conditions (jSDI), (j-^ljl and ()50|1 which correspond to the 
requirements (/cq, {X + iY) ci) = with f = I, Z ot ZZ. For n = 9 these are 
equivalent to 

02(38 + 70406 +70504 +0302 = (65a) 

50208 + 70406 —210004 — 7080-2 = (65b) 
20208 — 70406 -|-50802 = 0. (65c) 

These can be treated as a set of 3 linear equations in the 4 unknowns, 0203, 
O4O6, O6O4, O8O2 from which one finds that the group ()65p is equivalent to 

0203 = —aga2 = —ii^ (66a) 
3 

O4O6 = —aga^ = 1-1^ (66b) 

for some real parameter u. Note also that Re 0405 = Re 0203 = implies that 
any real solutions must have O4O6 = O2O8 = 0. However, all such solutions have 
been found above and none satisfy the additional requirements below. Hence, 
correcting all double-Z errors does require complex coefficients. 

Next we consider the conditions (jH^ . (jH^ and (j^H) which correspond to the 
requirements (/cq, (X — iY) Ci) = with f = I, Z or ZZ. These become 

ao08 + 280206 + 70|o4|^ +280502 +0309 = (67a) 
9ao03 + 1400206 + 70|o4p — 84o602 -70300 = (67b) 
9ao03 + 560206 — 70|o4p +50300 = (67c) 
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which can be rewritten as 



Reaotts + 28 Rea2a6+ 35|a4p =0 (68a) 
7 Re aoas + 28 Rea2a6— 35|a4p =0 (68b) 
Im aoOg + 14 Ima2a6 =0 (68c) 



since ()67ap and ()67b|) have the same real part while ()67b|) and ()67cp have the 
same imaginary part. Equations ()68a|) and ()68b|) are equivalent to 



35 

ReaoOs = — |a4p (69a) 
5 

Rea206 = — |fl4p- (69b) 
3 

To these conditions we need to add the requirements {Xcj^ i Ycj) = {cj, Zcj) = 
{Zcj, ZZcj) = which become. 

|ao|^ + 20|a2|^ + 14|a4|^-28|a6|^-7|a8|^ = (70a) 
9|ao|^ + 40|a2|^ - 14|a4|^ - 35|a8p = (70b) 

Im 00(22 + 21 Ima2a4 + 35 Im 0403 + 7Ima6a8 = 0. (70c) 

The last equation ()70cj) is precisely the condition Im {Xcj,iYcj) = and is 
obtained as the reduction of (j^J^ when n = 9 and conditions I and II hold. 

To recap, we have three groups of equations; namely, ()65|) from the conditions 
bfx + ibfY = 0, (jUTj) from the conditions bfx — ibfy = 0, and (fTUj) from the 
conditions dxy = diz = dz,zz = 0. In what follows, we will use the equivalent 
conditions in place of and or ()68c|) in place of ^7\i. 



6.5 Limits on correction of ZZ and one-bit errors 

To analyze the conditions obtained above, write = Xk + iyk- We can assume 
without loss of generality that Og = 1; then ()66ap implies that a2 = iv and 



implies Uq = — J;|a4p- We also have that ()66b|) implies Qq = Yj^o-ii from which 
we can conclude 2/6 = TWF' Equating the two expressions above for j/g yields 



= (71) 
9 X4 ^ ' 



which implies that X4 < 0. Under the assumption og = 1, ()70cp becomes 

7 



= z/ (xo — 21x4 + 15 ye 



V ( y |a4p + 211x41 + 15 + ^1^41') (72) 
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The last equation implies that either = or 04 = 0, either of which generates 
only a trivial solution. Thus there is no non-trivial solution to the seven equations 

dng), dnzj) and utib . 

By Theorem [7[ this implies that there is no 9-bit permutationally invariant 
code which can correct all one-bit errors as well as one type of double error. 

One might wonder if there is a 9-bit code which satisfies all the conditions 
above, except (jTUcJ- Such a code would still be of some interest. It would be able 
to correct all single and double Z errors, and detect all single X and Y errors. 
However, it would not be able to correct X^. and errors because it could not 
reliably distinguish between them. Unfortunately, even this is not possible. 

We return to the equations (j65j) and (|67p and observe that there are infinitely 
many solutions that can be expressed using one complex variable 04, or two 
real variables 0:4, j/4, in either case with the constraint Ima4 = X4 < 0. Let 
X = —Xi > and y = j/4. Then we have 

ao = f (1 - 2^y±) |a4p lao^ = + y^f 

6 X4 9 x"^ 

, .VSE |a4p , ,2 35(x2 + y2)2 

a2 = il^ = ±« ; 02 — 



3 V\^\ 9 



X 



a4 = X4 + iyi \ai\'^ = x^ + y'^ 



. v35 04 2 5 X + y 

ag = '^i—z 1= Ictel ^ 



T ^\^\ 7 X 

ag = 1 lagp = 1. 

Substituting into (fTOk ) and ()70b|) yields two equations in two unknowns which 
have no solution. Thus, there is no 9-bit code which satisfies all the desired 
equations except (jTOcj) . 

We also considered the possibility of dropping all Yk conditions to find a code 
which could correct all errors of the form single Xk, single Zk and double ZjZk. 
However, this is as restrictive as dropping only ()70c|) . 



7 Concluding Remarks 

Permutationally invariant codes which can correct all one-bit errors require a 
minimum of seven qubits. We have shown that there are two distinct 7-bit codes 
of this type. Although one might expect that 9-bit codes could also correct one 
class of double errors, a detailed analysis shows that this is not possible. Even 
a 9-bit code which could correct all one-bit errors of the form Xk and Zk and 
all two bit errors of the form ZjZk does not exist. If one modifies this to the 
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requirement that the code be able to correct all one bit errors of the form Xk and 
double errors of the form XjXi^ and ZjZ^, this can be done. However, it does 
not require 9-bits; it can be achieved using the simple 5-bit repetition code 
which can correct all double errors of the form XjX^ and YjYk as well. 

Permutationally invariant codes are highly degenerate, since all (2) exchange 
errors are equivalent to the identity. As discussed in Section l^TSl and illustrated by 
the 5-bit repetition code, this degeneracy can sometimes lead to enhanced ability 
to correct two-bit errors. However, there are also limitations on their ability to 
correct all two-bit errors of a given type as well as all one-bit errors, as shown 
by our analysis of of 9-bit codes. Although the reasons for this remain unclear, 
it may be that the "degeneracy enhancement" also gives hidden constraints, i.e., 
that one is implicitly trying to correct more two-bit errors than those from which 
the conditions were obtained. 

We have concentrated here on the construction of permutationally invariant 
codes. Actual implementation would require a number of additional considera- 
tions. For example, one would need a mechanism for initializing the computer 
in states corresponding to \cq) ® |co) . . . \cq). One could then obtain any state of 
the form \ck^) ® Icfcj) . . . |cfc„) with ki G {0, 1} by apphcation of for 
suitable choices of k. One also needs a mechanism for decoding, including a set 
of measurements which can distinguish between the different error subspaces, as 
well as a circuit for implementing the error correction process. 

Finally, one needs a set of gates for universal computation. As noted in 
Section 12.31 the logical X and Z operations are easily implemented as prod- 
ucts of their single-bit counterparts. The logical Y is given by the product 
i[®jXj\[®jZj\ = (-l)("-i)/2[^^.Y^.] ^]^gn ^ jg Moreover, the X,Y,Z gates 

all lie in the commutant of Sn- This implies that all gates needed for universal 
computation (including a non-trivial two-bit gate) also lie in the commutant since 
they can be written as linear combinations of these operations, the identity, and 
their products. (Alternatively, one could observe that all logical actions on code 
words must be permutationally symmetric and, hence, lie in the commutant of 
Sn-) However, this does not necessarily mean that all the desired gates can be 
implemented as products of a small set of one and two-bit gates. We leave the 
question of a practical implementation of a universal set of gates on code words 
for further investigation. 

We have only begun to explore the potential of non-Abelian stabilizer codes for 
quantum error correction; other examples should be studied. In addition to the 
issues identified above, there may be others which arise if one wants to combine 
non-Abelian stabilizers with other approaches to fault tolerant computation. 
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A Differences of one-bit errors 

In this section we will need some additional notation. Let Sr denote the binary 
n-tuple with components Ej = 6jr so that v + Er has components Vj + 6jr with 
addition mod 2. Let 1 be the binary n-tuple with all elements equal to 1. We will 
use s{v) = {j : vj ^ 0} to denote the support of of f = (t>i, . . . , w^,). 
It will be convenient to also introduce the vector 

\4(r,s)= {\^^^r)-\v\E,)) (73) 

wt(ti) = fc 
r,s ^ s{v) 

which is well-defined for k = 0, 1, ... (n — 2), and has the following properties 
when r s. 

{Vk{r,s),Wi) = for all r, s, A;, £ (74) 
{Vk{r,s),Vi{q,t)) = for m ^ £ V r, s, g, t (75) 

{V,{r-s),V,ir,s)) = 2 ^ ^) (76) 
(Vfc(r, s), Vk{r, t)) = ( ^ 7 ^ ) ^' ^' ^ all distinct. (77) 



k 

These are all straightforward, except (j77j) which follows from 



{Vk{r,s),Vk{r,t)) 



n — 3\ /n — 3 
k / ^ U - 1 



and the easily verified combinatoric identity (";,^) = ("fc^) + {^-l)- 

An important consequence of (fTUI) and (f77|) is that they imply that, for s ^ t, 

{Vk{r, s), \4(r, s)) = 2 {Vk{r, s), \4(r, t)). (78) 

This result plays an essential role in section 14.41 
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Our main results are that, for any code of the general form Q, 



n-l 



{Zr-Zs)\c^ = -2^afc Vfc„i(r, s) (79a) 

k=l 

n-2 

{Xr-X,)\co) = ^(afc-afe+2)Vfc(r,s) (79b) 

k=0 
n-2 

i{Yr-Ys)\co) = ^(afc + afc+2)Vfc(r,s) (79c) 

k=0 

with similar equations for |ci) and bk- Under the assumption that conditions (I) 
and (II) hold, we find the following variants useful 

{n-l)/2 

{Zr-Zs)\ci) = -2 ^ an^2m-iV2m{r,s) (80a) 

m=0 
(n-3)/2 

[{Xr-Xs)+t{Yr-Ys)]\Co) = 2 J2 (^2mV2m{r,s) (80b) 

m=0 
(n-l)/2 

[(X, - X,) - 2 {Yr - F,)] |co) = -2 ^ V2m-2{r, s). (80c) 

m=l 

To prove ()79ap and (IHUj) in the case of Z^, it suffices to observe that 

(Z, - Zs)Wu = I %y^_^^^^ ^) fo' ^ < ° < n - 1 

which is easily verified. 

Equations ()79b|) and (j79cp can be verified by some rather straightforward, but 
tedious, computations and combinatorics. One approach is to write out the effect 
of the errors Xr and Y^. Since these results are identical except for the signs of 
some terms, we introduce 

_ / +1 for X 
^xy = I _^ (82) 

and write the equations only for X^. with the understanding that these results 
hold for Yr with the sign changes indicated by ujxy 

XrWo = \er) 

XjWk = ^ \v + Er) + uJxY ^ \u) for 1 < k < n — 1 

wt(ti) = k wt(ji) = fc — 1 

r ^ s(t;) r ^ s{u) 

XrWn = UJxY\l + £r)- 
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For distinct r, s, we want to determine the effect of the differences Xr — X^, 
Yj. — Yg on the Wk, and for this purpose, the following expression, which we write 
only for 2<A;<r;, — 2, is useful. 

XrWk = ^ \V + Er) + ^ \V + Eg + Er) 

wt{v) = k wt(i') = fc — 1 

r,s ^ s{v) f,s^ s(v) 

+ ^XY ^ \u) + UJXY ^ \U + E,) (83) 

wt{u) = fc — 1 wt{M) = fe — 2 

r,s (fi s(u) T^s ^ s{u) 

Then (|H^ implies 

{Xr-Xs)Wo = \Er)-\Es)=Vo{r,s) 

{Xr-Xs)Wi = J2{\E,+Er)-\E,+Es))=Vi{r,s) 

{Xr-Xs)Wk = Vk{r,s)-u;xYVk-2{r,s) 

{Xr-Xs)Wn-l = UJxY "^ill + £j + Sr) -\l+Ej + £s)) = -i^XYVn-3{r,s) 
{Xr-Xs)Wn = UJxYi\l + er) -\l+Es)) = -0OxYVn-2{r,s). 



To analyze double phase errors, first observe that when (/q, ZZcj) = (with 
/ = X,Y,Z), the analogous inner products involving (n — l)-dimensional rep- 
resentations will be zero if and only if ((/i — ft)ci, ZrZsCj) = 0. By consider- 
ing the action of the transposition (It), one can show that this holds whenever 
{1, t} = {r, s} or {1, t} n {r, s} = 0. Hence, it suffices to consider r = 1 and t s, 
in which case one can use (I79ap and ZiZg = / — ZxiZx — Z^ to conclude that 

n— 1 

ZiZ,|co) = |co) + 2^afcZi\4_i(l,s). (84) 
fc=i 

We will also need the formula 

(....(M,,..(M,)^(:-)-("-v^("^^) (-) 

which follows from (|7j) with X = n — 

We now let f = Z. Using (|79a|) again and (|85p with k = 2m — 1, one finds 
{{Zi — Zt)co, ZiZgCo) = if and only if 

m=l ^ 
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Then it follows from ^(IZl) = T(~-i)ln:2) il) that §^ is equivalent to (l56b)l 
minus Ti^—n times equation ()56a|) . Thus, the "block" conditions for double phase 
errors do not add additional constraints when conditions (I) and (II) hold. 

The cases / = X ± zy, and z = 0, j = 1, can be dealt with similarly, but are 
not needed here. We note only that, unlike the case f = Z, they do generate 
additional constraints. 

B Decomposition into irreducibles 

In view of the repeated use of decompositions into irreducible subspaces, we 
explicitly write out some of them. Recall that C^" = Q)^=QWk', each Wk is the 
eigenspace of the operator Sz = \ J2k = with eigenvalue ^{n — 2k). Each 
Wk can be further decomposed into a direct sum of subspaces which transform as 
irreducible representations of Sn- In fact, the action of Sn on Wk or Wn-k is its 
action on sets of size in {1, . . . , n}. For < < L|J , Wk is known to decompose 
into a sum of irreducible subspaces indexed by the partitions [n — j, j] for j = 
0, . . . ,k, each appearing once. Physicists may recognize that this is equivalent to 
a decomposition into simultaneous eigenspaces of and the total spin^ operator 
S2 = Si + S^ + si with eigenvalue s(s + 1) for s = f , f - 1, . . . i|n - 2k\. 

To facilitate use in counting arguments, as in Section EH we adopt the con- 
vention of labeling an irreducible subspace (in part) by its dimension. Thus, 
denotes an irreducible subspace of Wk with dimension d. 

For n = 5 



Wo 


= Ul 






Wi 








W2 


= uli 






W3 


= ul^ 






W4 


= uU 






W5 


= l^l 







^In this one paragraph, we use the famihar Sx,Sy,Sz, rather than the equivalent 
^X, ^Y, ^Z, to denote the components of spin, and trust that context suffices to distinguish 
them from the symmetric group denoted Sn, which is a very different entity. 
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For n = 7 



yvQ 


— LAq 










= Uli 








W2 


= Uli 








m 


= l^u 








W4 










W5 


= ul^ 








m 


= uli 








W7 


= u] 









where U2 denotes the irreducible representation associated with the partition 
[5, 2] and Ul^ a second, distinct, 14- dimensional irreducible representation asso- 
ciated with [4, 3]. In terms of spin, has s = | and Ul'^ has s = |. 

For n = 9 



Wo 


= 1^0 










Wi 


= Uli 










W2 


= U\i 










W3 


= Ul^ 


dUl^ 








W4 


= Uli 


^Uli 




D 6 




W5 


= uli 


dUli 








We 


= l^l^ 










W7 












Ws 


= ul^ 










Wg 


= 1^9 











C Complex coefficients 

If the ttk are not real, then one must modify the analysis in Section 0] accordingly, 
and require both real and imaginary parts of the resulting equations to be zero. 
We again use the classification of error conditions described at the end of Section|Sl 
We omit the details and summarize the results. 
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a) Setting the real parts of (jHUjl - ljHHI) to zero yields the pair of equations 

(n-l)/2 



= ^4^(^)«i±i +2 k(^^Re{akan-k+i) (87) 

n + 1 / n \ (n-3)/2 
= +2 5^ (n- A;)r j Re(afca„_fc_i). (88) 



The imaginary parts of both ()3()|1 and (jH^ are always zero and do not place 
any additional restrictions on Ofc. Setting the imaginary parts of ()31|) and 
(jnni) to zero yields the conditions 



(n-l)/2 



^ Im[afca„_fc+i] /c(n-2/c + l) f j 
fc=i ^ ^ 

= ^ lm[akan-k^i]{n- k){n~2k~l) i y 

b) The condition ()38|1 from the off diagonal terms in D becomes 

= J2\^k\\n-2k) (j). 

1 — n ^ / 



(89) 
(90) 



(91) 



(This is sufficient to ensure that djz = and Djz = 0, as well as that the 
real part of (ji^ is zero.) To ensure that the imaginary part of {Xci,iYci) 
is zero we must also require 



n — 2 



k=l 



= ^Im[afc+iafc_i].(^_ J. (92) 



(This also ensures that the imaginary part of PHj) is zero so that dxY = 
and DxY = 0.) 

As before, we analyze the "block" conditions only under the assumption that 
conditions (I) and (II) hold. As for real coefficients, these conditions do not yield 
new requirements. 

c) Setting the imaginary parts of (j40p and (j4H) to zero yields conditions equiv- 
alent to (IHni) and (TOjl . 

d) The expression in (jl^ gives two conditions. The first is equivalent to 
and the second to (j^^ with k = 2m + 1. 
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